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Abstract: Interest rate market models, like the LIBOR market model, 
have the advantage that the basic model quantities are directly observ¬ 
able in financial markets. Inflation market models extend this approach 
to inflation markets, where zero-coupon and year-on-year inflation-indexed 
swaps are the basic observable products. For inflation market models con¬ 
sidered so far closed formulas exist for only one type of swap, but not for 
both. The model in this paper uses affine processes in such a way that 
prices for both types of swaps can be calculated explicitly. Furthermore 
call and put options on both types of swap rates can be calculated using 
one-dimensional Fourier inversion formulas. Using the derived formulas we 
present an example calibration to market data. 


Arbitrage-free inflation models have first been rigorously introduced in Jarrow and Yildirim 
[10]. Since then several inflation models have been proposed. Similar to interest rate models 
one can distinguish between short rate models and market models. While short rate models in 
the spirit of Jarrow and Yildirim [10] aim at modeling the unobservable continuous nominal 
and real short rate, market models use discrete observable rates as the basis for modeling 
(see Belgrade et al. [l], Mercurio [16]). These observable rates are the basis of liquidly traded 
inflation swaps, zero coupon inflation-indexed swaps and year-on-year inflation-indexed swaps. 
While there have been several extensions of these models (e.g. Mercurio and Moreni [17, 18]) 
all of these models suffer the problem that there exist analytical formulas for only one type 
of swap, but not both. The model in this paper leads to closed formulas for both types. 

Based on the ideas in Keller-Ressel et al. [14] one can use affine processes to describe analyt¬ 
ically highly tractable models. Affine processes are Markov process, where the characteristic 
function is of exponentially affine form, i.e. 

E [e“^*|As] = 

The class of affine processes contains a large number of processes, e.g. every Levy process is 
affine. Using the longest-dated nominal zero coupon bond as a numeraire we can model the 
normalized bond prices as “exponential martingales” with respect to an affine process X. In 
this type of models we are not only able to price both types of inflation swaps, but can also 
derive semianalytical formulas for calls and puts on the underlying inflation rates, another 
liquidly traded inflation derivative. 

The structure of this paper is as follows. The first part describes inflation markets and 
typical traded derivatives. Afterwards we outline the setup of an inflation market model. 
The second part introduces the affine inffation market model and derives pricing formulas for 
the introduced inflation derivatives. The third part provides a concrete model specification 
including the calibration to actual market data. In the appendix we collect the properties of 
affine processes needed in this paper. Furthermore we specify the affine processes used in the 
numerical section. 
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1. Inflation markets 


Denote the time t price of the nominal zero coupon bond with maturity T by P{t,T) and 
consider an inflation index with time t value I{t). Typically inflation indexes are so-called 
consumer price indexes (CPI). To shorten notation we will use the term CPI synonymous 
for inflation index, nevertheless the reader can think of an arbitrary inflation index. The 
basic mathematical instruments in inflation-linked markets are so called inflation-linked zero 
coupon bonds (corresponding to zero coupon bonds for nominal interest rate markets). An 
inflation-linked zero coupon bond with maturity T is a bond paying I{T) at time T. Denote 
its price by PiLB{t,T). 

In actual markets governments issue inflation-linked coupon bonds. Such a bond pays a 
hxed coupon on the variable basis I{Tk)/I{T q) at some fixed number of predetermined dates 
Tk <T (typically annually), where Tq is the time of issue. Additionally to the coupons such 
a bond redeems at maturity T with value max{/(r)//(To), 1}. Such a bond can therefore 
be described as a combination of inflation-linked zero coupon bonds plus an included option 
with payoff (1 — /(T)//(To))+. In general these bonds are issued with maturities of several 
years and inflation is positive. In this case the included option has little influence on the total 
price, which is why it is market practice to mostly ignore it. In particular if one ignores these 
options it is possible to strip inflation-linked zero coupon bond prices out of actually traded 
inflation-linked coupon bonds by the same methods used for nominal quantities. 

Consider the quantity 

Pn(t.T) := ( 1 ) 

which is called the price of a real zero-coupon bond. Note that this is not the price^ of a traded 
asset, but a theoretical quantity. The usage of the term price is motivated by the fact that this 
quantity can be viewed as the price of a zero coupon bond in a hctitious economy, where ev¬ 
erything is measured in terms of the inflation index /(t)^. Given real zero-coupon bond prices 
continuously compounded real interest rates are dehned by R{t, T) := — In {PR{t, T))/{T — t). 
Accordingly one can define real counterparts to other nominal quantities such as forward in¬ 
terest rates or the short rate. This quantities are sometimes used as a starting point for 
inflation option pricing models (see e.g. Jarrow and Yildirim [10], Mercurio [16]). 

Next to inflation-linked bond markets there exist several liquidly traded inflation-linked 
derivatives. First consider the forward price of the inflation index (forward CPI), i.e. the at 
time t fixed value I{t, T), which at time T can be exchanged against I{T) without additional 
costs. Since PiLsit^ T) is the current price of /(T), the time t forward CPI for maturity T is 


x(t,r) 


PlLB{t,T) 

P{t,T) 


( 2 ) 


In a zero coupon inflation-indexed swap (ZCIIS) two parties exchange the realized inflation 
against a fixed amount (1 -|- ZCIIS are mostly traded for full year maturities M. 


^Here we mean price in terms of money which has to be paid at a transaction. In fact it is essentially this 
number that is quoted on trading screens. However, in case such a bond is traded, the cash-flow is then 
the qnoted number multiplied by the according index ratio. 

^One could interpret I{t) as a numeraire, but one has to be careful not to use this as a mathematical 
numeraire, since 7(f) is not actually traded. 
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For T = t + M the value of such a payer swap can be expressed as 




( 3 ) 


The rate K, for which equation (3) is zero is then called the ZCIIS rate ZCIIS{t; M). These 
ZCIIS rates are quoted in the market for several full-year maturities. 

Remark: Note that ZCIIS rates and inflation-linked bonds are closely related via (2) and (3). 
In reality this relationship is not observed. This is partly due to different creditworthiness of 
counterparties in bond and swap markets. A more detailed analysis of this difference can be 
found in Fleckenstein et al. [7]. For model calibration one has to choose one market, usually 
the swap market. 

Next to ZCIIS there is a second important type of swap in inflation markets, the year-on- 
year inflation-indexed swaps (YYIIS). These swaps exchange the annualized inflation against 
a fixed rate K, i.e. consider an annually spaced tenor structure Tk = t + k, k = 0,..., M. The 
netted payment of a payer YYIIS at time is (I(Tfc)//(Tfc_i) — 1) — K. Hence the inflation 
leg consists of payoffs of the form 


1 

T-S 



Denote the forward value of such a payoff, the annualized forward inflation rate, by Fi{t, S, T). 
Then the value of a payer YYIIS with maturity M and strike K can be expressed as 


M 

Y^P{t,Tk)iFiit,Tk_i,Tk) - K). 

k=l 


The YYIIS rate YYIIS{t\M) is the rate K such that the corresponding YYIIS has zero 
value. Note that given YYIIS rates for all annual maturities one can calculate annual forward 
inflation rates Fj{t,Tk-i,Tk) and vice versa. 

Toward CPIs respectively forward inflation rates Fj(t,S,T) are the mathematical 

quantities underlying the market-traded ZCIIS, respectively YYIIS. Inflation market mod¬ 
els aim at modeling these quantities. In existing inflation market models either X(t, T) 
or Fj{t,S,T) can be expressed by analytical formulas, but not both. Consider a market, 
where price processes are assumed to be semimartingales on a hltered probability space 
(D, (77), P). Fix a T-forward measure i.e. a probability measure equivalent to P 

such that asset prices normalized with the numeraire price P{t,T) are Q^-martingales. Then 


Fi{t,S,T) 


[iiT)\:Ft 

1 

T-S 


m 

I{S) 




Calculating the expectations of the inflation index, as well as the fraction of the inflation 
index at two different times proves difficult. In the model of this paper both are exponentially 
affine in the underlying driving stochastic process. For an affine process (see appendix) such 
expectations can be calculated and we are able to give semianalytical formulas for many 
standard options like caps and floors of forward inflation rates. 
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1.1. The inflation market model 


We now introduce the general setup of an (inflation) market model. Consider a tenor structure 
0 < Ti < ■ ■ ■ < T]\f =: T and a market consisting of zero coupon bonds with maturities 
Tfc and prices P{t,Tk). The price processes (P(t, Tfc))o<t<Tfc are assumed to be positive 
semimartingales on a filtered probability space (fl, {^t)o<t<T, IP), which satisfy P{Tk,Tk) = 
1 almost surely. If there exists an equivalent probability measure such that the normalized 
bond price processes P{-,Tk)/P{-,T) are martingales^, the market is arbitrage-free. This 
setup describes the class of interest rate market models like the classical LIBOR market 
model (Brace et al. [2]) and its extensions. 

To extend this setup to inflation markets consider an inflation index I, where we assume 
w.l.o.g. that 1(0) = 1. Assume there exist inflation-linked zero-coupon bonds with the same 
maturities'^ Ti,... ,Tn and price processes PiLB{t,Tk)o<t<Tj,, all of which are positive semi¬ 
martingales®. If there exists an equivalent probability measure such that all normalized 
price processes 

( Pit,Tk) \ ( PlLB{t,Tk) \ 

\ J 0<t<Tk \ P{t^P) ) 0<t<Tk 

are Q^-martingales, the extended market model is arbitrage-free. For given define the 
Tfc-forward measures by 


dQ^'' _ 1 P{0,T) 

dQ^ ~ P{Tk,T) P{0,Tky 


Under the forward interest rate 


F^{t) 


1 ( P{t,Tk-i) \ 

Ak\ P{t,Tk) )' 


the earlier introduced forward CPI 


Afc := Tk — Tfc-i, 


P{t,Tk) 


PlLB{t,Tk) 

p{t,n) 


and for j < k the forward inflation rates Fj{t,Tj,Ty given by 


1 + (Tfc - Tj)Fj{t, Tj,Tk) 




are all martingales. Modeling (some of) these martingales is the starting point of inflation 
market models in the literature (see e.g. Mercurio [16]). In contrast we start by modeling the 
normalized bond prices in (4) and derive the above quantities thereof. 


^One can extend bond price processes to [0,r] by setting P{t,Tk) := for t > Lc, so that 

P{-,Tk)/P{-,T) is a martingale on [0,T] if and only if it is a martingale on [OjTfc]. Economically this 
can be interpreted as immediately investing the payoff of a zero coupon bond into the longest-running zero 
coupon bond. 

'^The assumption that for each zero coupon maturity there is a ILB with the same maturity is used only for 
notional convenience. 

®The inflation index is only described through the bond prices Pile- I.e. the distribution of I{t) is only 
given at times T*, where it coincides with the distribution of PiLB{Tk,Tk)- 
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2. The AfFine inflation market model 


Let {Xt)Q<t<T with Xq = x be an analytic affine process with state space M>q x M"', m > 0, 
n > 0 on the probability space (n, A, {Xt)o<t<T, Q^) and define for = 1,..., 


P{t,T) ’ 

PlLB{t,Tk) __ 
P{t,T) ’ 


Uk e (M^ox{or)n V, 
Vk £ n V, 


( 6 ) 


where 

M“ := [e“-^^ | Pt] = exp (.^T-i(«) + V’T-t(n) • X*) , u£V, (7) 

with V defined in (19). The processes are Q^-martingales by the definition of an affine 
process. Hence this model is arbitrage-free. Note that in (6) the parts of corresponding 
to real-valued components of X are chosen to be zero. For a decreasing sequence ui > • • • > 
UN P 0 one then has > M“'', so that forward interest rates F^{t) are guaranteed 

to be nonnegative for all k. Contrary to interest rates® inflation rates are not required 
to be nonnegative, which is why we do not restrict in (6). The values of and 
should be calibrated to fit the initial term structures, i.e. = P{0,Tk)/P{0,T) and 

Mq*" = P/2.b( 0, Tfc)/P(0, T). By Lemma 2 it follows that parameters fitting a current term 
structure with nonnegative forward interest rates can always be chosen to be decreasing. For 
multidimensional affine processes such sequences are far from unique. Concrete specifications 
how to choose Uk and Vk will be presented in section 3. 

The big advantage of this setup is that (5) in this case reads 


1 




exp 


{4>T-n{uk) + ipT-niuk) ■ ^Tfe) 


( 8 ) 


which is exponentially affine in X. In particular it is easy to check (see Keller-Ressel et al. 
[14]) that for 0 < s < r and il)T-r{uk) + w £V 


I =exp(^(/)r_s(V’T-r(tifc) +w) - 4>r-si'lpT-7 


exp 


^-sii’T-riUk) +W) - V’r-s('0T-r('Ufc))) ' 77*^. 


(9) 


Hence the moment generating function of X is also known' under different measures QP. 
Together with the exponential affine form of basic quantities this is the reason why this model 
is analytically highly tractable. For example the forward rates satisfy 


{l + AkF\t)) 


t — A{t,u,,_^,Uf,)+B{t,uu-i,Uk)-Xt 


^Although interest rates are currently negative in certain countries, interest rates are still bounded below by 
the costs of physically keeping money. We can incorporate bounds different from 0 by setting ~ 

CkMp. 

'"This also shows that X is a time-inhomogeneous affine process under . 
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with 


( 10 ) 


A{t,v,u) := (t)T-t{v) - (j)T-t{u), 

B{t, V, u) := ipT-t{v) - ipT-tiu). 

Hence the -extended moment generating function of In (l -|- AkF^{t)) can be calculated 
explicitly using (9). The price of a caplet then follows using a Fourier-inversion formula (see 
Keller-Ressel et al. [14]). Swaptions can also be dealt with (Keller-Ressel et al. [14], Grbac 
et al. [9]) and so the most common interest rates derivatives can be calculated efficiently. We 
can use similar methods for inflation derivatives. 


2.1. Forward CPI and CPI options 

As mentioned before, the main advantage of this model is that for several important quantities 
the moment generating function is known under all forward measures Start by looking 
at the forward CPI 




PlLB{t,Tk) 

P{t,Tk) 


PlLBit,Tk) P{t,T) ^ ^ A{t,v,,,Uk)+B{t,Vk,UkyXt 

P{t,T) P{t,Tk) 


( 11 ) 


with A and B defined in (10). Hence the forward CPI is of exponential affine form and 
therefore the -moment generating function of its logartihm can be calculated using (9). 
In particular, setting A’} := A{Tk,Vk,Uk), Bj := B{Tk,Vk,Uk) and using I{Tk) = I{Tk,Tk) 
one has 


M 




(z) := 


[iiTkr\Ps] = 


explzA'} + zB’} ■ Xt,)\P. 


= exp(^zA^j -b (pTk-sii’T-Tkiuk) + zBj) - 4>Tk_silpT-Tk 


exp(^{'lljTi^_silpT-niuk) + zBj) - IpTk-sii’T-niuk))) ■ 

= exp(^z4>T-TkiVk) + (1 - z)(pT-Tk{Uk)^ 
exp(^(l)Ti^_s{z'tpT-Tkivk) + (1 - z)'ll)T-n{uk))^ 
exp(^Tk-s[zi)T-Tk{vk) + (1 - z)ipT-n{uk)) ' Xs'j/M^'^. 


Here the last equality follows using (20). Note that this function is well defined and analytic 
in z if zi^T-nivk) + (1 - z)'ilJT-Tk{uk) £ int(V). 

Civen the moment generating function of ln(/(Tfc)) CPI calls and puts can be calculated 
using the following well-known Fourier inversion formula (see e.g. Eberlein et al. [5]). If 
R G (l,oo) such that Aix\jr{R) < oo, then 

( /r-(w+^) \ 

EHe-- KU^] ^ - I + (12) 

Thus the price of a forward CPI call with maturity Tk and payoff {I{Tk) — K)^ is 

KP(t Tu) r°° ( nT \ 

CPICalK*. n.K) = R, (» + R) du. 
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where i? > 1 is chosen to satisfy RipT-Tki'^Jk) + (1 “ R)'^T-Tj^iuk) £ int(V). 

Remark: In section 1 it is mentioned that ILBs usually come with an included option guar¬ 
anteeing a redemption of at least the original nominal amount. For an ILB issued at S and 
with maturity this translates into an option {l — I{Tk)/I{S))^ which corresponds to 1/I{S) 
CPI puts with strike I{S). 


2.2. Forward Inflation and inflation caplets and floorlets 

Typically inflation market models are not able to handle both forward CPI and forward 
inflation products analytically. With the approach presented here forward inflation rates are 
of a similar form as forward CPIs. The annualized inflation Fj{Tk,Tk-j,Tk) satishes 


1 + (Tfc - Tk-j)Fi{Tk, Tk-j,Tk) 


m-j) ■ ■ 


(13) 


The following Lemma gives the moment generating functions for random variables of this 
type. 

Lemma 1: Let s < r < t < T and 


'll)T-t{Uk) +w eV and 'ijjt-ri'lpT-tiUk) + w) - 1pT-r{Uk) + uev. 


Then 


jgQ k I ^ eX.p(^[lpr-s{'lpt-r{'lpT-t{uk) + w) + u) - IpT-siUk)] ' 

exp(^(j)t-ri'‘pT-tiuk) +W) + 4>r-silpt-ri'4’T-tiuk) + w) + u) - (pt-sii’T-t 
Proof. Using the tower property and applying (9) twice it follows 


U-JCr-l-LO-Xt 1 77 

<3 *^5 


^u-Xr jr^ 


exp^(ft_ 

ri'lpT-tiUk) +W) - ft-ri'f’T-tiUk))'^ 


exp\^{'tpt-r{'f’T-tiuk) +w) - ift-ri'lpT-tiuk)) + u) ■ Xrj \ Fs 
= exp(^(l)t-r{'lpT-t{Uk) +W) - ft-ri'f’T-tiUk))^ 

exp(^(j)r-s{'f’t-r{'4’T-t{uk) + w) + u) - fr-sif^T-riuk))^ 
exp(^{'lfr-si'llJt-rilpT-t{uk) + w) + u) - ifr-sidpT-riuk))) ' X 
= exp(^(l)t-ri'4>T-tiUk) +W) + (pr-si'lpt-ri'f’T-tiUk) + w) + u) - ft-si'f’T-t 
r-silpt-rilpT-tiuk) + w) + u) - tpT-siuk)) ' . 


exp 


(14) 


Here we used the semiflow property (20) to simplify the expression. 


□ 
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By Lemma 1 the -moment generating function of defined in (13) is 


M 


Yk\T: 


{z) = 


exp(zAf + zB’^j • Xt, - zA]-^ - zB]-^ • | F, 

= exp(^zAf - zA)~^ + (pT^_Tk_j{i;T-n{uk) + zB'})'^ 
ew[(Pn-j-sii^n-n-j{'>pT-niuk) + zBf) - zb)~^) - (pT^^sitPr-T, 
expf(V’Tfe_,-s(V'Tfc-Tfc_,.(V'r-Tfc(^^fc) + zB’P) - zBj~^) - IpT-siUk)) ■ xX 


which is well-defined if 

l^lpT-n(Uk) + zB’} ,'lpTk-Tk-ii'lljT-TkiUk) + zB’}) - 'ipT-n.j (Uk) “ zBj~^] C V. (15) 

The forward inflation rate is then given by 1 -|- (T^. — Tk-j)Fj{t,Tk-j,Tk) = So 

for ipT-n-j{vk) - G V it is 

l + {Tk- Tk-j)Fi{t,Tk-j,Tk) = exp(^{'iljTk-j-t{'pT-Tk_j{vk) - B^^) - 'ipT-t{uk)) ■ Xt'^ 
exp(^(pT-Tk-j{uk-j) + (pTk-j-ti'tpT-Tk-j (Vk) - Bj •’) -h (pT-t 
Furthermore the payoff of an inflation caplet with strike K is 

(Tk-Tk-j)iFjiTk,Tk.„Tk)-K)+=(^X^^ 


- K 


where LC = 1 -I- (T^ — Tk-j)K. With the Fourier inversion formula (12) one can calculate the 
price of an inflation caplet. In particular, we have for i? > 1 

provided that ipT-niuk) + RB'} G int(V) and ipTk-n-jii’T-Tkiuk) + RBj) - ipT-n-jiuk) - 

rbX e mt(V). 


2.3. Correlation 

So far we considered the pricing of typical market traded options. Another important aspect 
is the correlation structure. The relevant quantities 


ln(l + AfcF„^(f)) = In 
In (X'^ (f)) = In 


M, 


Uk-l 


Ml 


= A{t, Uk-i,Uk) + B{t, Uk-i,Uk)Xt, 


A{t, Vj,Uj) + B{t, Vj,Uj)Xt. 


(16) 


ln(l -F AkFj{t,Tk-j,Tk)) = const -h {ipTk-j-MT-nivk) - b’} - ipT-tiuk)) ■ Xt. 
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are all affine transformation of Xt and the correlation for two such terms is 


Cor[^i + Bt • Xt, At + Bf Xt] 


Var[^t • Xt, Bt • Xt] 

^ Var[5t • Xt] ^ Var[5t • Xt] 


For independent components of Xt this simplifies to® 

Eti bib ; Var[X|] 

\/Etl(B;)" Var[X-] Var[X-] 

Hence correlations strongly depend on i?(t, Ufc_i, n^) = 'ijjT-tiuk-i)—'4’T-tiuk) and B{t, Vj, Uj) = 
'tpT-tivj) — 'ipT-t{uj), respectively the structure of the Vk and Uk- The exact correlation de¬ 
pends on the used measure (e.g. but choosing Vk and Uk cleverly, one can guarantee 

that the correlation structure, i.e. the correlation signs, stay the same. Concrete specifications 
for meaningful correlation structures will be given in the next section. Similar observations 
can also be made for instantaneous correlations of the corresponding quantities in the case of 
continuous affine processes (see Grbac et al. [9] for the general idea). 


3. Implementation example 

We design the structure of the affine inflation market model in such a way that the calibration 
can be separated into the calibration to nominal market data and the calibration to inflation 
market data afterwards. The method used to calibrate to nominal market data is based on 
the ideas in Grbac et al. [9]. There they fit a multiple curve affine LIBOR market model 
by using a common driving process X*^ plus additional driving processes X^,... ,X^, all of 
which are independent, affine and nonnegative. For calibration they use caplets with full-year 
maturities and underlying forwards of tenors less than a year. Using one individual driving 
process for each year one can then use an iterative procedure to calibrate to market data. 
Their approach can also be used in this setup. 

In particular consider a semiannual tenor structure = k/2,k = 1,...,X, N even, and 
a driving affine process consisting of M -|- 1 = N/2 + 1 components X*^,X^,... ,X^, all of 
which are independent analytic affine processes with functions (jf and V'®, f = 0,..., M. Then 
by (23) 


M 


Mu) = 

1=0 

Mu) = {'iptM)MlM),---,Miu^)), 


where u^,i = 0,... ,M denotes the corresponding component of u G To describe an 

affine inflation market model we also have to specify the vectors Uk- The structure of the 


®Up to some technical conditions the variance of a one-dimensional affine process is 


Var[W] 


d'^u 


u=Q 
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Ui 

U2 

U3 

Ui 


UN-2 

UN-l 

UN 


X^ 

X^ 



XM 

j^M+1 

j^M+2 


j^2M 

Ul 

Ul 

Ms 


UN-l 

0 

0 


0 

U2 

U2 

U3 


UN-l 

0 

0 


0 

U3 

0 

U3 


UN-l 

0 

0 


0 

Ui 

0 

Ui 


UN-l 

0 

0 


0 

UN-2 

0 

0 


UN-l 

0 

0 


0 

UN-l 

0 

0 


UN-l 

0 

0 


0 

UN 

0 

0 


UN 

0 

0 


0 


Table 1: Description of the parameter structure Uk- Each row corresponds to one vector with 
the column names denoting the process the position in the vector corresponds to. 


vectors should be so that the following points are satisfied. 

• Forward interest rates are nonnegative. This is the case if 0 < < Uk-i- 

• The model matches the initial interest rate term structure. This basically fixes one 
component of each vector Uk- 

• Calibration to market data is possible by using an iterative procedure. 

• The model has a meaningful correlation structure. 

This can be achieved by choosing the vectors in the following way. They depend on 2N 
real parameters 

Ul > • • • > UN >0, Ui > • • • > UN > 0. 

For 1 < j < M set (compare table 1, ignoring the zero columns in the rightmost side for now) 

M 

Uk = UkfP+ Uke\^'^ + ^ U2i-ie^ 

«=rii+i 


where e°,e^,...,e^ denote the base vectors (1,0,..., 0),..., (0,..., 0,1) of Note 

that with this choice the vectors (uk) are decreasing. Given ui,...,Un and the processes 
X^, X ^,..., X^ by Lemma 2 the parameters ui,..., un are determined by fitting the cur¬ 
rent term structure. I.e. we require that 


P{0,Tk) 

p(o,r) 






r - vr^il 

M 

^U2l-\Xlp 




-ni+i 



so that the parameters ui,... ,un can be calculated using backwards iteration. Furthermore 
the semiannual forward interest rates with full-year maturities only depend on U 2 k-i,U 2 k 

and the processes X^,X^,... ,X^. Hence if X^ and ui,... ,un are already specified, one 
can fit X^ to caplets on the forward interest rate and then go backwards to iteratively 
fit the processes to caplets on forward interest rates Hence if X^ and ui,... ,un 

are fixed, all the remaining parameters can be calibrated to the yield curve and caplet prices. 
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Figure 1: Yield curve and ZCIIS curve from September, 29th 2011 


As stated in Grbac et al. [9] and confirmed by our numerical tests the concrete choice of 
and Uk (within a meaningful range) has no qualitative impact on the resulting calibration 
quality. Henceforth is fixed as a CIR process (specifications of this process can be found 
in the appendix, see equation (24)). So far we have said nothing about resulting correlations. 
This is where the choice of the parameters ui,... ,Un comes in. The relevant functions for 
correlations of forward interest rates are 

B{t, U2k-l,U2k) = {i>T-t{u2k-l) - V’T-t(^‘2fc), 0, . . . , 0, '4)T-t{u2k-l) “ ''/'T-t(«2fc), 0, . . . , 0). 

These functions are “orthogonal” except for the first component®. Hence by equation (17) 
the correlation structure mainly depends on the sequence (ufc). Since A® is nonnegative the 
function tptiu) is increasing in u (see Keller-Ressel et al. [14]). With (uk) being a decreasing 
sequence this results in nonnegative correlations. Setting Uk = u for all k results in zero 
correlation^®. The magnitude of correlations depends on how much of the variance of forward 
bond prices is explained by UkX^. Here we choose the factors Uk so that E*® = 

P{0,Tk)/P{0,T). The idea behind this choice is that approximately half of the variance 
should be explained by the common factor. Alternative if one has additional information on 
correlations (e.g. through market data), this could be incorporated in Uk- 
For calibration we used market data from September, 29th 2011. The yield curve is boot¬ 
strapped from LIBOR and swap rates and is displayed in figure 1. The used caplet implied 
volatilities for 6-month forward interest rates with strikes 1% to 6% and maturities from 1 
to 10 years are bootstrapped from cap data. The resulting implied volatilities can be found 
in figure 2. We hxed the time horizon T = 10 and chose the parameters of the common CIR 
process as A = 0.026, 6 = 0.65, r] = 0.5, x = 3.45. The M individual driving processes are cho¬ 
sen to be CIR processes with added jumps (see appendix, equation (25)). Their parameters 


®Although only stated for even forward rates for notional simplicity this is true for all forward rates. 
^'^Negative correlations in this setup are only possible if the sequence (uk) is not decreasing which means that 
forward interest rates can become negative. 
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Figure 2: Market and model caplet implied volatilities of 6-month forward interest 
rates with strikes 1% to 6% and maturities from 1 to 10 years. Market volatilities 
are in a transparent blue, while model volatilities are displayed in red. 


and the parameters Uk are calibrated with the mentioned recursive method. In each step the 
parameters are chosen so that the mean squared errors of implied volatilities are minimized. 
In contradiction with Grbac et al. [9] we were not able to produce a similar calibration quality 
as described in their paper^^. The resulting calibration can be found in figure 2. Especially 
for long dated caplet volatilities the pronounced skew could not be reproduced. Nevertheless 
the model provides a reasonable fit of the caplet volatility surface, especially since the focus 
is on inflation derivatives. 

The next step is to extend the calibration to inflation markets. Additionally to the M + 1 
driving processes used for modeling interest rates we use another M independent analytic 
affine processes (one for each year) driving inflation-related quantities. Since the individual 
processes are independent by (23) the setup without inflation can be embedded by setting 
the additional components of Uk equal to zero (see table 1). So from now on we assume that 
the processes X ^,...,and the vectors ui,..., are already fixed. 

The choice of the inflation parameters Vk focuses on the same aspects aspects as the choice 
of Uk without the restriction that inflation rates have to be nonnegative. We again assume 
that the vectors Vk are determined by 2N parameters'^ ■ ,vn, vi,. .. ,vn- In particular, 
we choose (see also table 2) 


Several requests for clarification with the authors have resulted in the answer that their results are currently 
not in a state to be shared. 

^^We actually only require N parameters, since the odd rows in table 2 do not contribute for the considered 
annual inflation rates. For notional simplicity we nevertheless consider 2N parameters. 
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Table 2: Description of the inflation parameter structure Each row corresponds to one vector 
with the column names denoting the process the position in the vector corresponds to. Note 
that for annual inflation rate option pricing the vectors Vj with j odd do not matter. 


M 

Vk = Vke°+ Uke^^'\ + ^ U2i-ie''+ Vke^~^^^'\. 

^=m+i 

Choosing the nominal components = u^. for i = 1,...,M has the advantage that for¬ 
ward CPIs and therefore also forward inflation rates do not depend on the nominal pro¬ 
cesses X^,... Together with the choice of Vk with respect to the inflation processes 

.. .X"^^ this implies that the forward CPI X{t,T 2 k) depends only on X^ and 
In particular the function B{t, V 2 k, U 2 k) corresponding to the forward CPI X(t, T 2 k) defined in 
(11) is 

B{t, V2k, U2k) = (V’T-tfefc) - 'tpT-tiu2k),0, . . . , 0, 0, . . . , 0). 

From this it follows that for different forward CPIs these functions are “orthogonal” except for 
the first component. Furthermore except for the first component they are also “orthogonal” to 
the functions B{t, Uj-i, Uj) relevant for forward interest rates. Hence the correlation structure 
mainly depends on Uk and Vk- Since is monotonically increasing, one should choose Vk > Uk 
if the corresponding forward CPI should be positively correlated with nominal interest rates 
or Vk < Uk if it should be negatively correlated^^. Also note that two annual forward CPIs are 
positively correlated if sgn (ufc — Uk) = sgn (hj — uj). This therefore gives us some criteria 
how to determine the parameters Uk from correlation assumptions and from now on we assume 
that the parameters vi,...,vn are given. Assuming also a fixed process X^^^ we would 
like to determine Vk from the current term structure PiLB{0,Tk)/P{0,T). We differentiate 
between two cases. First consider a nonnegative affine process In this case by Lemma 

2 there exists a unique Vk, so that Mq*" = PiLB{0,Tk)/P{0,T). For Vk ~ Uk this typically 
results in Vk > 0. In this case zero coupon inflation for [t, Tk] is always positive. To avoid this 
alternatively consider an affine process, which takes negative and positive values. In this case 
Mq is not necessarily increasing in Vk- However, by Lemma 2 it is still a convex function in 
Vk- This means that there are at most two possible choices for Vk, in which case we need to 

decided to introduce correlation for the inflation part only via the common factor process. We could 
also have changed the parameters Uk in order to introduce correlation. In this case even more complicated 
correlation patterns could be created. 
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Figure 3: Linear interpolated annual forward inflation rates Fj(0, T2(k-l),T2k) 
(black) and its approximation I{t, T 2 k)/F{t, T 2 (^k-i)) ~ 1 (dashed red) for matu¬ 
rities from 1 to 10 years. 


pick one. For the results in this paper we only used results where one choice was smaller than 
zero and the other one larger than zero, which we then picked. To determine the processes 
,..., notice that the annual forward inflation rate F/(t, T 2 k- 2 ,T 2 k) depends only 
the processes , X^ and Fj{t, Tq, T 2 ) depends only on X^, By starting with 

inflation options on Fj{T 2 ,Tq,T 2 ) one can calibrate the parameters of X^^^. Then one can 
iteratively calibrate the parameters of X^~^^ using inflation options on Fi{T 2 k-,T 2 [k-i)-,'F 2 k)- 

For the calibration example we used ZCIIS rates from 1 to 10 years (see figure 1) and 
inflation options for strikes ranging from —2% to 6%, the prices of which are displayed in 
figure 4. We chose Vk = 11^(1 + ck) with c ~ 0.08, so that ck is between 1 and 1.15. This 
means that correlation are positive as usually observed. For the processes ... ,X‘^^ 

we used Ornstein-Uhlenbeck processes with added jumps (see appendix, equation (26)). The 
parameters of these processes and the sequence (vk) are then calibrated as described. In each 
step the mean squared errors of option prices are minimized. The resulting fit is displayed in 
figure 4 and this shows that the calibration is very accurate. 

We would also like to display the fit in terms of (shifted lognormal) implied volatilities. Typ¬ 
ically one only has quotes on ZCIIS rates, which do not directly translate to forward inflation 
rates. However, annual forward inflation rates can be approximated by Fj(t,Tk-j,Tk) ~ 
I{t, Tk)/F{t, Tk-j) — 1 (see figure 3). It is market practice to use this as forward value in the 
shifted Black formula to calculate approximate market volatilities. The resulting fit in terms 
of shifted implied volatilities is then displayed in figure 5. This shows that implied volatilities 
are closely reproduced by the model across all maturities and smiles and that the model is 
very capable in fitting the observed market data. 

Conclusion 

We have introduced a highly tractable inflation market model, where we are able to derive 
analytical formulas for both types of inflation-indexed swaps. Furthermore inflation caps 
and floors as well as CPI caps and floors can be calculated with a one-dimensional Fourier 
inversion formula. Hence prices for liquidly traded inflation derivatives can be calculated 
quickly and accurately. Additionally the proposed model is able to price classical interest 
rate derivatives like caps and floors. Using these formulas we are able to calibrate the model 
to market data. The calibration example shows that the model can be calibrated to inflation 
market data very accurately. 
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Figure 4: Market and model caplet/floorlet prices in basis points for annual forward 
inflation with strikes —2% to 6% and maturities from 1 to 10 years. Market 
prices are in a transparent blue, while model prices are displayed in red. For 
strikes between —2% and 1% prices are quoted for floorlets, for strikes between 
2% and 6% prices are quoted for caplets. Market prices are bootstrapped from 
corresponding cap/floor data. 



Figure 5: Market and model implied volatilities for annual forward inflation with 
strikes —2% to 6% and maturities from 1 to 10 years. Market volatilities are in 
a transparent blue, while model volatilities are displayed in red. 
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A. AfFine processes 


Let X = {Xt)o<t<T be a homogeneous Markov process with values in D = M>q x M"' on 
a measurable space (O, M) with filtration {Xt)o<t<T, with regards to which X is adapted. 
Denote by and E* [•] the corresponding probability and expectation when Xq = x. X is 
said to be an affine process, if its characteristic function has the form 

jgx _ exp^(/)i(u) + V’t('w) • u€ iM'^,x G D, 

where (/> : [0,T] x IM'^ C and V' : [0,T] x C'^ and iM'^ = {u £ : Re(u) = 0}. By 

homogeneity and the Markov property the conditional characteristic function satisfies 

E^ jT^] = eyip(4>t-siu) + ipt-siu) ■ Xs^ (18) 

Accordingly affine processes can also be defined for inhomogeneous Markov processes (see 
Filipovic [6]). In this case the affine property reads 

E® [e“'^‘I X's] = exp(^(j)s^tiu)+ 'ips,t{u) ■ Xs'^, u£iW^,x£D, 

with : iM'^ —)■ C and 'ips,t : —)• for 0 < s < t. 

X is called an analytic affine process, if X is stochastically continuous and the interior of 
the set^^ 


V := 


r 


\u£&: sup E"^ 

^Re(w)-Xs 

1 o<s<r 

I- 


< oo 


Mx £ D 


(19) 


contains 0. In this case the functions (p and ip have continuous extensions to V, which are 
analytic in the interior, such that (18) holds for all u G V (see Keller-Ressel [11]). 

The class of affine processes includes Brownian motion and more generally all Levy processes. 
Since Levy processes have stationary independent increments, it follows that 'iptiu) = u, while 
(pt{u) = tK{u), where k is the cumulant generating function of the Levy process. Ornstein- 
Uhlenbeck processes are further important examples of affine processes. The affine processes 
used in this work are described at the end of this section. 

The standard reference for affine processes is DufRe et al. [4]. There they give a characteri¬ 
zation of affine processes, where p and ip are specified as solutions of a system of differential 
equations. To motivate this consider an affine process X. By the tower property for condi¬ 
tional expectations it holds for all x G D 




= E"= [E"^ [e^^^+pXs]] = E^ 


^4>t{u)+'>pt{u)-Xs 


can be described as the (convex) set, where the extended moment generating function of Xt is defined for 
all times t and all starting values x. By Lemma 4.2 in Keller-Ressel and Mayerhofer [12] the set V is in 
fact equal to the seemingly smaller set |m € C'* : G int(_D) : IgUsi^i VTl ^ qqI 
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Using equation (18) it follows that (j) and satisfy the so-called semi-flow equations 


(l)t+s{u) = Mu) + (j)s{Mu)), (po{u) = 0 , 

Msiu) = ijjsitptiu)), Mu) = U. 


( 20 ) 


For a stochastically continuous affine process X it was shown in Keller-Ressel et al. [13] that 
the functions 


F{u) ■■= ^Mu] 


t=o+ 


R{u) ■■= ^Mu) 


i=0+ 


exist^®. Rewriting (20) in terms of difference quotients and letting s —?■ 0 we get that (j) and 
if: satisfy generalized Ricatti equations 



= F{Mu)), 

Mu) = 0, 


= R{Tpt{u)), 

Mu) = u. 


( 21 ) 


The functions F and R have a specific form of Levy-Khintchine type as first described in 
Duffie et al. [4]. There it is also shown that for every F and R of this form (21) has a unique 
solution. Specifying the functions F and R is an alternative way to specify an affine process. 
Keller-Ressel and Mayerhofer [12] give conditions on F and R under which a solution of (21) 
defines an analytic affine process. Note that in order to evaluate 4> ijj one would like to 
have closed form solutions to the system (21), which in general is not the case. 

Coupling independent affine processes is very tractable. For two independent affine processes 
X and Y and all starting values x, y one obtains 




^(ux,UY)iXt,Yt} 




( 22 ) 


Hence {X,Y) is an affine process with 

4^’^\ux,uy) = (l>f{ux) + <t>l {uy), 

'iI^\^'^\ux-,uy) = {M{ux)AJ {uy)). 

This fact together with the following Lemma is used in section 3. 

Lemma 2: Let X he an analytic affine processes comprised of m + n independent affine 
processes, where the first m are nonnegative. For {ui, ..., Uk,v, Uk+i, ..., Un) G int(V)n]R™'“’“"' 
define the function 

Then f^ is monotonically increasing if k < m and f^ is convex for all k. 

Proof. For each x £ D the term inside the expectation is convex in v and monotonically 
increasing in u if fc < m. This then also holds after taking the expectation. □ 

The last part of this section describes the affine processes used in this paper. One classical 



^®This was also shown for affine processes with general state spaces in Keller-Ressel et al. [15] and Cuchiero 
and Teichmann [.'Ij. 
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example is the CIR process, which is the unique solution of 


dXt = -X{Xt -9)dt + 27? dWt, Xo = X. 


(24) 


For this process the functions (j) and tp are defined for Re(u) < ^(1 — e 


Mu) = - 


Mu) = 


xe 

2r?2 


In 1 - 


2rf 




e 




-w 


u 


The CIR process almost surely stays nonnegative. It is strictly positive if ^ > t?^. One can 
add jumps to this process by adding the differential of a compound Poisson process Lt to the 
dynamics of X. 


dXt = -X{Xt -e)dt + 2rj^t dITt + dL*, 


Xn = X. 


(25) 


If Lt has exponentially distributed jumps with expectation values ^ arriving at rate A/3 the 
functions (p and pj are (see Grbac and Papapantoleon [8]) 


Mu) = - 


Mu) = 


xe 

2rf 


In ( 1 — 


2rf 


{l-e-Mu - 


A/3 

A — 2rfa 


In 


a — u 


a 


— 7x (e + (1 — e 


e 


l-^(l-e-^‘)n 


where 




1 


Re(u) < min ■( ;^(1 — e ^*) ^,a(e + (1 — e ^^)^2'rf‘a 


-1 


,a 


Since Lt has only positive jumps this process also stays nonnegative. As a third example 
consider the real-valued affine process defined by 


dXt = -X{Xt -e)dt +a dWt + dLt, 


Xo = X, 


(26) 


where Lt is a compound Poisson process with positive jumps with mean ^ arriving at rate 
A/3"*“ and negative jumps with mean ^ arriving at rate A/3“. The functions (p and ip in this 
case read (see Miiller and Waldenberger [19]) 


Mu) = 

+ 


aM 


4A 

/3+-/3 


(1 - e-M + - e-M + In 


[a+ - e“^* 


u){a -|-e ^*u) 


(a+ — u) {a + u) 


In 


(a"*" — e ^^u){a + u) 
(a+ — u){a~ + e~^^u) J ’ 


Mu) = e 

for —a~ < Re(7i) < a~^. 
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